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We investigate ®Li(n, n')^Li* — >■ d + a reactions by using the continuum-discretized coupled-channels 
method with the complex Jeukenne-Lejeune-Mahaux effective nucleon-nucleon interaction. In this study, the 
^Li nucleus is described as a d + a cluster model. The calculated elastic cross sections for incident energies 
between 7.47 and 24.0 MeV are good agreement with experimental data. Furthermore, we show the neutron 
spectra to ®Li breakup states measured at selected angular points and incident energies can be also reproduced 
systematically. 

PACS numbers: 24.10.Eq, 25.40.Fq 
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I. INTRODUCTION 

The ^Li nucleus is known to have a well developed d + 
a cluster stracture and the binding energy is very small 1.47 
MeV from the d-a breakup threshold. Because of those fea- 
tures, a breakup process of ^Li into d and a is one of the 
significant mechanisms of the reaction, and the accurate de- 
scription is required. As one of the most reliable methods for 
treating breakup processes, the continuum-discretized cou- 
pled channels (CDCC) method | lil] has been proposed and suc- 
cessfully applied to analyses of three-body breakup systems, 
in which the projectile breaks up into two constituents, such 
as ^Li into d and a. Thus, the CDCC method is expected to 
be indispensable to analyze ®Li breakup reactions. 

For the CDCC analyses of ^Li elastic and inelastic scatter- 
ing on various targets, Sakuragi et al. iQ] have suggested to 
reproduce the absolute values of a lot of experimental data of 
cross sections very well. They have also resolved the anoma- 
lous renormalization problem [i3J| in folding model potentials 
of ^Li by showing large dynamical polarization potentials due 
to breakup processes. Furthermore the CDCC method has 
been also applied to the sequential Coulomb/nuclear breakup 
via the resonance state of ^Li, in which both Coulomb and 
nuclear breakup processes were taken into account consis- 
tently ji!]. The Coulomb breakup is of a great interest as the 
time-reversed reaction of the radiative capture reaction, d-\- a 
— > ^Li H- 7, which is one of a key reaction in the nucleosyn- 
thesis in the early Universe or during stellar evolution ^ . So, 
it is important that the method could describe not only the nu- 
clear breakup process but also the Coulomb breakup one on an 
equal theoretical footing. In spite of their efforts, the CDCC 
method has not been confirmed yet its applicability to breakup 
to continuum states of ^Li, especially breakup spectra of ^Li. 
This is a basic interest of the theoretical framework to describe 
reaction mechanisms relevant to breakup processes. 
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The n H- ^Li reactions are important not only from the ba- 
sic research interest as shown above but also from the ap- 
plication point of view. Lithium isotopes will be used as a 
tritium-breeding material in d-t fusion reactors, and accurate 
nuclear data are required for n- and p-induced reactions. In- 
deed, IAEA is organizing a research coordination meeting to 
prepare nuclear data libraries for advanced fusion devices, 
FENDL-3, and the maximum incident energy is set at 150 
MeV to comply fully with the requirements for the IFMIF de- 
sign 16:;]. Lithium isotopes are ones of important materials in 
these libraries. As a matter of fact, the neutron interaction 
with lithium and associated breakup reactions are indispens- 
able in determining the neutron energy spectra in blankets of 
fusion reactors. Thus the tritium breeding ratio, nuclear heat- 
ing distributions and radiation damage of structural materials 
are affected by the n-n Li reactions significantly. 

In spite of the importance of the n H- ^Li reaction as de- 
scribed above, experimental data leading to ^Li continuum 
breakup processes are extremely rare for the neutron energy 
region above 20 MeV |0-|9tl- Furthermore, the statistical 
model, used often in evaluation of nuclear data for medium 
to heavy nuclei, cannot be applied to the ^Li(n, n') reactions. 
As mentioned in Refs. ifiol [Till , the mechanisms leading to 
three- or four-body final states, n -\- d -\- a and n -\- n -\- p -\- 
a, are important. Therefore more reliable theoretical calcula- 
tions for the cross sections are highly desirable. 

In this paper, cross sections for ^Li(n, n')^Li* d -\- a 
reactions are evaluated by using the microscopic coupled- 
channels method |0, [l2l, in which we adopt microscopic 
wave functions of ^Li with the d -\- a model and the complex 
Jeukenne-Lej£une-Mahaux effective nucleon-nucleon (JLM) 
interaction [13] between n and ^Li. The calculated elastic, 
inelastic cross sections and continuum neutron spectra to ®Li 
breakup states show good agreement with the experimental 
data at low neutron incident energies, where breakup effects 
of ^Li into d and a are dominant. Furthermore we discuss ap- 
plicability of the approach used here to the present system for 
the neutron energy region up to 150 MeV. 

This paper is organized as follows. In Sec. |II] we describe 
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the CDCC method including calculations of ^Li wave func- 
tions and coupling potentials with the JLM interaction. In 
Sec. Hn] and Sec. IIVI we present the calculated results and 
discuss the applicability, respectively. Finally, we give a sum- 
mary in Sec.rv] 



II. FORMULATION 

The n + ^Li scattering system is described as a n + d + a 
three-body model, and the Schrodinger equation is written by 



Kr- 



6 Li 



E 



^itR) = 0, (1) 



where R and ^ represent a coordinate between n and the 
center-of-mass of ^Li and a set of internal coordinates of ^Li, 
respectively. Kn is a kinetic energy operator associated with 
R, i/aLi is the internal Hamiltonian of ^Li, and an effective 
interaction between jth nucleon in ^Li and n is represented 
by Vjo- 

The total wave function with the total angular momentum 
J and its projection M on z-axis, '^jm, is expanded in terms 
of the orthonormal set of eigenstates of -ffsLi- Here, micro- 
scopic wave functions of ^Li based on the d-a cluster model 
are written as 



(2) 



where r is the relative coordinate between d and a. The d- 
a relative wave function with the angular momentum £, total 
spin / and its projection m on z-axis, ?/'/™, can be described 
as 
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with the spin of deuteron, j]d = 1. The internal ground state 
wave functions of d and a, (p{d) and (p{a), are assumed as 
inert core. In the present calculation, we use the same in- 
ternal model Hamiltonian of the d-a system, Hej^^, adopted 
in Ref. |2tl, which are fixed to reproduce the experimental 
data of the ground state (-1.47 MeV) and 3+-resonance state 
(0.7 1 MeV) as well as low-energy part of d-a scattering phase 
shifts. The ground state is assumed to be a pure S state 
{£ = 0), and the calculated wave functions well reproduce 
not only the elastic electron scattering form factor but also the 
inelastic one to the 3+-resonance state |2|]. In this framework, 
we take into account the totally antisymetrized effect by as- 
suming the orthogonal condition model ifTill . 

In the CDCC method, the relative wave functions of ^Li 
for the ground and continuum states are described as a finite 
number of discretized states denoted by 

4'irir) = ^iAr) [z%{nr) ^ {i = O-N), (4) 

whose energies ejf, are given by 



(V'5"(r)|i/eLi|V'5"(r))r. 



As the discretization approach, we adopt the pseudo-state 
method ifTsl - flTll here. The advantage of the pseudostate 
method is that if there are resonances in its excitation spec- 
trum, we can describe discretized continuum states with a rea- 
sonable number of the basis functions, without distinguishing 
the resonance states from non-resonant continuous states as 
mentioned in Ref. lfl6ll . As well known, ^Li has three reso- 
nances in £ ~ 2. Therefore the pseudostate method is very 
useful for analyses of ^Li breakup reactions. 

In the pseudostate method, we diagonalize Hei^^ in a space 
spanned by a finite number of L^-type basis functions for 
r, {<Pji}, and discretized wave functions for the radial part, 
{(t>l(} are obtained by 



Hiir) = £A},^,,(r). 



(6) 



As the basis functions, the complex-range Gaussian basis 
functions lfl6l [isll are adopted, and we include £ ^ and 2 
for the relative angular momentum of d-a: ^ S represents the 
5-wave state with / = 1, and ^D, ^D, and "^D correspond 
to the D-wave state with 1 — 1,2, and 3, respectively. The 
number of states N is decided as involving all open channels 
for each incident energy. 

After the discretization and the truncation of d-a contin- 
uum, 5* jA/ is reduced to an approximate one. 
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with 



yjM = [[^%{^r) ®Vd]j«> l''YL{nR)]jM^{a)ip{d), (8) 

where 7o = (0, 0, 1, L, J) and 7 = (i, L, J) represent 
the elastic channel and breakup channels, respectively. The 
expansion-coefficient Xj in Eq. (Q represents the relative mo- 
tion between n and ^Li, and L is the orbital angular momen- 
tum regarding R. The relative momentum is determined 
by the conservation of the total energy. 



E = P^/2fi + ei,. 



(9) 



(5) 



with II the reduced mass of the n + ^Li system. Multiplying 
Eq. ([Til by ('")3^jm] * from left, one can obtain a set of cou- 
pled differential equations for x-y, called the CDCC equation. 
Solving the CDCC equation under the appropriate asymptotic 
boundary condition, we can obtain the elastic and discrete 
breakup S'-matrix elements. Details of the formalism of the 
CDCC method are shown in Ref. 

For the calculation of diagonal and coupling potentials in 
the CDCC equation, we use the complex Jeukennd-Lejeune- 
Mahaux effective nucleon-nucleon (JLM) interaction lfl3ll 
based on a single folding model. The JLM interaction can 
be easily applied to coupled-channels calculations of nucleon- 
nucleus systems, such as Refs. lfl2l[l9tl . This interaction has 
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FIG. 1 : Angular distribution of the elastic differential cross section of n + Li scattering for incident energies between 7.47 and 24.0 MeV. The 
solid and dashed lines correspond to the result with and without couplings to breakup states of ^Li, respectively. Experimental data are taken 
from Refs. Q-QSl- The data are subsequently shifted downward by a factor of 10^^-10"^ from 8.96 MeV to 12.0 MeV on the left panel, and 
10^^-10"" from 13.9 MeV to 24.0 MeV on the right panel, respectively. 



energy and density dependence, and forms a Gaussian with 
both real and imaginary parts; 

Vjo{Rjo]p,E) = A,l/(p, £;) cxp (-i?joAfl) 

+ iXn^Wip, E) exp i-Rjo/tj) ,(10) 

where Rjq is a coordinate of a nucleon in ^Li and n. The pa- 
rameters, tfl, tj, and Ay, are taken as the same ones used in the 
original paper iflSll . = ti = 1.2 and ~ 1.0. Meanwhile 
the normalization for the imaginary part, A^,, is optimized to 
reproduce the elastic cross sections, because the strength cor- 
responding to the loss of flux depends on the model space 
considered in the calculation. Details for strengths of V{p, E) 
and W{p, E) are shown in Ref. iflsd . 

Using the JLM interaction, the diagonal and coupling po- 
tentials, Ky-y', are obtained by 

^77'(^) = J Py','{s,^R)vjo{E,p,Rjo)dsdnR, Ci-i) 

where transition densities, p^^', and averaged matter density, 
p, of ^Li between 7 and 7' are defined by 
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FIG. 2: Neutron inelastic scattering angular distribution to the 3 - 
resonance state 0.71 MeV above the d-a threshold. The experimental 
data are taken from Ref. |[l^. The data for 14.1 (18.0) MeV are 
shifted by a factor of 0.1 (0.01). 



III. RESULTS 



p{s) = - j {pyj{s,ilji) + py'y'{s,ilii)} dftsd^lji, (13) Figure [T] shows the differential elastic cross sections of 



respectively. Here Sj is a coordinate of jth particle in ^Li 
relative to the center-of-mass of ^Li. 



+ °Li for incident energies between 7.47 and 24.0 MeV. One 
sees that the results of the CDCC calculation represented by 
the solid lines are in good agreement with the experimental 
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FIG. 3: The neutron spectra calculated by the CDCC method with the JLM interaction comparing with measured data at selected angular 
points in the laboratory system. The experimental data are taken from Ref. fioll . 



data. The dashed lines represent results of a single channel 
calculation, in which couplings to breakup states are omitted. 
It is found that breakup effects shown by the difference be- 
tween the dashed and solid lines are significant to reproduce 
the angular distributions of the elastic scattering. For all in- 
cident energies, we take A^, = 0.1 to reproduce the data. The 



main component of the flux loss from the elastic channel is 
due to breakups of ^Li to d and a, which can be taken into 
account in the CDCC calculation directly. The other effects 
of the flux loss, such as the excitation of a and d breakup, 
are not so large because of low incident energies. Therefore 
the strength of A^, becomes very small. Here, it should be 
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noted that the single channel calculation cannot reproduce the 
experimental data if any values of are taken. 

Figure ID shows the angular distributions to the first excited 
3+ state of ^Li for E,, =1 1.5, 14.1, and 18.0 MeV. The theo- 
retical cross sections are calculated by integrating the breakup 
cross section to 3+-continuum for the resonance energy re- 
gion. One sees that the CDCC calculation can also reproduce 
the inelastic cross section. 

In Fig. [3] the calculated neutron spectra are compared with 
the experimental data at selected angles in the laboratory 
frame and incident energies. Components of ^5, ^D, ^D, 
and are represented by the dash-dotted, dashed, dotted and 
thin solid lines, respectively, and these results are broadened 
by considering the finite resolution of the experimental ap- 
paratus fioll . Three peaks in the experimental data represent 
the elastic, inelastic to the S'*" -resonance, and 2+ -resonance 
components, respectively, from higher neutron energies. The 
CDCC calculation gives a good agreement with experimental 
data at the higher neutron energy region. On the other hand in 
the low neutron energy region, which corresponds to high ex- 
ited states of ^Li, the calculated cross section underestimates 
the experimental data for all energies and angles. The result 
indicates that experimental data contain contribution from the 
(n, 2n) reaction corresponding to a four-body breakup reac- 
tion ^Li(n, nnp)a, as mentioned in Ref. ifioll . In the present 
calculation, the four-body breakup effect is not taken into ac- 
count directly, and the effect is treated as a absorption effect 
on the elastic cross section. 
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FIG. 4: Neutron energy dependence of the total cross section of n + 
®Li scattering. The dotted line represents the experimental data ll20h . 
The open circles, squares, and triangles corresponds to results with 
Xw = 0.1, 1.0, and 1.5, respectively. 



higher incident energies (> 30 MeV) there is no experimental 
data of the elastic cross section, which are needed to optimize 
the normalization factor of the imaginary part. One of the op- 
timization methods is to fit the neutron total cross section data 
of ^Li, which have been measured up to 560 MeV ||2Q|] . Fig- 
ure]?] shows the measured total cross section comparing with 
the calculated ones with A„, =0.1, 1.0, and 1.5 represented 
by open circles, squares, and triangles, respectively. One sees 
that large values for A^, are required with respect to increas- 
ing the neutron energy. Although the energy dependence for 
A„, is rationalized because absorption effects become strong 
as increasing the energy, the optimized value at 150 MeV is 
much large about A„,=1.5. This normalization problem would 
be resolved by analyzing p + ^Li elastic cross sections, which 
have been measured at higher incident energies, within the 
present framework. Furthermore for higher incident neutron 
energies four-body breakup effects of ®Li(n, nnp)a are not 
negligible. Therefore we require the four-body CDCC calcu- 
lation, which has been successful for analyses of ^He breakup 
reactions ll2llj25ll . in the higher energy region. 



V. SUMMARY 

We analyze ^Li(n, n')^Li* d + a reactions by us- 
ing the continuum-discretized coupled-channels method with 
the complex Jeukenne-Lejeune-Mahaux effective nucleon- 
nucleon interaction. In the present analysis, it is found that the 
elastic cross sections for incident energies between 7.47 and 
24.0 MeV can be reproduced by the present analysis with one 
normalization parameter for the imaginary part of the JLM in- 
teraction [Xu, = 0.1), and breakup effects on the elastic cross 
section are significant. Furthermore the calculated inelastic 
cross section to the 3+ -resonance state and neutron spectra 
are also good agreement with the experimental data system- 
atically. Thus, the CDCC method with the JLM interaction 
is expected to be indispensable for the data evaluation of the 
^Li(n, n') reactions, and the advantage is to obtain not only 
elastic and inelastic cross sections but also neutron spectra 
within the same framework. 

For higher neutron incident energies, we have also dis- 
cussed the applicability of the present microscopic approach. 
It is found that the required normalization factor is much 
large, X^j = 1.5 for 150 MeV, from the analyses of the to- 
tal cross section. In order to investigate the normalization 
problem we should analyze p + ^Li elastic scattering, which 
have been measured at higher incident energies. Furthermore 
the four-body CDCC calculation is required to describe the 
four-body breakup process of ^Li into n + p + a, which be- 
comes significant for the higher energy. These results will be 
reported in a forthcoming paper 



IV. DISCUSSION 

Theoretically the CDCC calculation with the JLM interac- 
tion can estimate cross sections of the n + ^Li scattering for 
neutron energies up to 160 MeV, which is the maximum appli- 
cable energy suggested in the original paper lfl3l . However for 
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